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■ ^ ', Abstract 

CZ2 ! 

We present first-principles total-energy calculations of (001) surfaces of SrTiOs. 
■y ■ Both SrO-terminated and Ti02-terminated surfaces are considered, and the results 
C \ are compared with previous calculations for BaTi03 surfaces. The major differences 
are in the details of the relaxed surface structures. Our calculations argue against 
^ ■ the existence of a large ferroelectric relaxation in the surface layer, as had been 
^ , previously proposed. We do find some indications of a weak surface ferroelectric 
instability, but so weak as to be easily destroyed by thermal fluctuations except 
Q ■ perhaps at quite low temperatures. We also compute surface relaxation energies 
and surface electronic band structures, obtaining results that are generally similar 
to those for BaTi03. 
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^ 1 Introduction 

CN ■ 
O ' 

oo : 

^ ' The cubic perovskites are an important class of materials that are of particular interest because 
"ctf ■ of the variety of structural phase transitions that they display [1]. The structural instabilities 
may be of ferroelectric (FE) character, as for BaTiC>3, or of antiferrodistortive character (in- 
volving rotation of oxygen octahedra), as for SrTi0 3 . Specifically, SrTi0 3 adopts a paraelectric 
£h ! simple cubic perovskite structure above T C =105K, but transforms to a tetragonal antiferrodis- 
^ ; tortive structure below T c ; and while it shows strong signs of FE fluctuations at very low 
■ temperature, it evidently remains paraelectric down to zero temperature. 

• i-H , 

^ \ The (001) and (111) surfaces of cubic perovskites have been most investigated experimentally 
& ' [2]. For a II-IV perovskite such as SrTiC>3, there are two possible nonpolar (001) surface ter- 
minations: a SrO-terminated surface (type-I), and a TiCVterminated surface (type-II). On the 
other hand, the (111) surfaces are polar, and therefore presumably much less stable. Here we 
focus on the (001) surfaces of SrTi0 3 . 

Interest in the surface properties of SrTi0 3 arises because of the catalytic properties of these 
surfaces [3], and because of the common use of SrTi0 3 as a substrate for epitaxial growth of 
high-T c superconductors (such as YBa 2 Cu 3 7 ) [4] and other oxides. 
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Theoretical studies of these surface have been numerous. Wolfram and coauthors [5], using a 
linear combination of atomic orbitals approach, predicted mid-gap surface states for SrTiC>3, 
in disagreement with experimental investigations [6,7]. Tsukada et al. [8] employed the DV Xa 
cluster method to study SrTi0 3 surfaces, finding no mid-gap surface states. However cluster 
methods are not very suitable for high-accuracy calculations of relaxations and electronic states 
on infinite surfaces, underlining the need for the application of more accurate, self-consistent 
techniques. First-principles density-functional calculations have been very successful in the 
study of bulk perovskites [9-12], and more recently there have been similar calculations for 
perovskite surfaces [13-16]. In particular, Kimura et al. [16] applied the same method as in this 
paper, the plane- wave ultrasoft-pseudopotential method [17], to the study of (001) surface of 
SrTiOs, with and without oxygen vacancies at the surface. The main difference between that 
report and the present work is that we fully relax the atomic coordinates in the slab. Also, 
we analyze the possible existence of a FE surface layer for the SrO surface [18], as has been 
suggested previously. 

On the experimental side, the study of these surfaces is complicated by the presence of surface 
defects [19], making it difficult to verify the surface stoichiometry. The results also tend to 
depend upon the surface treatment [20]. Low-energy electron diffraction (LEED) and reflection 
high-energy electron diffraction (RHEED) studies of SrTi0 3 surfaces have been reported [20,21]. 
These show no reconstruction of the surface layer. In contrast, photoelectron spectroscopy [22] 
and scanning tunneling microscopy [23] observations show different reconstructions for the 
reduced surfaces. The absence of in-gap surface states has also been shown [6]. 



2 Theoretical approach 

As in previous work on BaTiC>3 surfaces [13], we employ here a self-consistent pseudopoten- 
tial technique in which the valence electron wavefunctions are obtained by minimizing the 
Kohn-Sham total-energy functional using a conjugate- gradient technique [24]. The exchange- 
correlation potential is treated within the local-density approximation (LDA) in the Ceperley- 
Alder form [25]. Vanderbilt ultrasoft pseudopotentials [17] are used to avoid norm-conservation 
constraints, thus permitting the use of a small plane-wave cutoff of 25 Ry, in spite of the fact 
of that we are dealing with first-row and transition metal atoms. Such a cutoff has previously 
been shown to be adequate in the bulk [24]. The forces on each ion were relaxed to less than 
0.02 eV/A using a modified Broyden scheme [26] . 

Our calculations are carried out in a periodic slab geometry. For the type-I (SrO terminated) 
surface, the slab contains 17 atoms (4 SrO layers and 3 Ti0 2 layers). Similarly, the type-II (Ti0 2 
terminated) slab contains 18 atoms (4 Ti0 2 layers and 3 SrO layers). In both cases, the slabs 
were thus three lattice constants thick; the vacuum region was two lattice constants thick. The z- 
axis is taken as normal to the surface. The calculations were done with a (6,6,2) Monkhorst-Pack 
mesh [27]. This k-point set produced results of very good accuracy. The structure was set up 
using our theoretical lattice constant of 3.86 A, which is about 1% smaller than the experimental 
one; this underestimation is typical of LDA calculations. For further details concerning our 
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Table 1 

Equilibrium atomic displacements (relative to ideal positions) for the SrO- and Ti02-terminated sur- 
faces, when no symmetry-breaking distortions are allowed. Units are at the theoretical lattice constant 
(a = 3.86 A). 



Layer 


SrO Surface 


8z 


Ti0 2 Surface 


Sz 


1 


Sr 


-0.057 


Ti 


-0.034 




Oi 


0.001 


On 


-0.016 


2 


Ti 


0.012 


Sr 


0.025 




On 


0. 


Oi 


-0.005 


3 


Sr 


-0.012 


Ti 


-0.007 




Oi 


-0.001 


0„ 


-0.005 



method and its accuracy, we direct the reader to our earlier paper [13]. 

Below 105K, SrTi0 3 transforms to a tetragonal antiferrodistortive (AFD) state in which the 
oxygen octahedra rotate about a (001) axis in opposite directions in alternate unit cells [28]. 
However, these rotations are typically small; even at T = they amount to only ~ 3°. Since we 
are mainly interested in comparing with room-temperature experiments, and since the static 
AFD distortions have already disappeared at 105K (well below room temperature), we have 
not included them in our calculations. On the other hand, we do want to consider the possible 
presence of a FE surface layer, as suggest in Ref. [18]. 

We thus choose symmetries as follows. Most of the calculations were carried out using a "full set" 
of symmetries, consisting of M x , M y , and M z mirror symmetries (the latter being a reflection 
through the center plane of the slab), as well as 1 x 1 translational symmetries parallel to the 
surface. The full set of symmetries prevents the occurrence of FE as well as AFD distortions. 
Then, some further calculations were carried out using a "reduced set" of symmetries, identical 
to the full set except that M x symmetry is allowed to be broken. The reduced set of symmetries, 
while still suppressing the AFD instability, allows the surface to develop a FE distortion along 
x if it should turn out to be energetically favorable. 



3 Surface relaxations 

We begin by presenting the relaxed structure for each of the two surface terminations, obtained 
by starting from the ideal structure and then relaxing the atomic positions while preserving 
the full set of symmetries. The relaxed geometries are summarized in Table 1. (Coordinates are 
only listed for atoms in the top half of the slab, z > 0; the others are determined by M z mirror 
symmetry). By symmetry, there are no forces along x or y. To set notation, oxygen atoms Oi 
are the ones lying in SrO planes of the slab, while On refers to the oxygen atoms lying in Ti0 2 
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Table 2 

Comparison of theoretical and experimental structural parameters for symmetry-preserving surface 
relaxations. Ad\2 and Ad 2 3 are respectively the changes in interlayer spacing for the first and second 
pair of layers, while s measures the outward displacement of the surface oxygens relative to the first- 
layer metal atoms. All quantities are in A. 







Q 






DL \J~ LUI lllllldloU. 












n 


99 




n i n 


Theory, Ref. [29] 


0. 


.14 


-0.47 




Expt., Ref. [21] 


0. 


.16 


-0.19 


0.08 


Expt., Ref. [20] 


0. 


.16 


0.10 


0.05 


Ti02-terminated 










Theory, present 


0. 


.07 


-0.27 


0.12 


Theory, Ref. [29] 


0. 


.04 


-0.39 




Expt., Ref. [21] 


0. 


.08 


0.04 


-0.04 


Expt., Ref. [20] 


0. 


.10 


0.07 


0.05 



planes of the slab. The layer numbering in Table 1 is from the outermost layer inwards. 

From Table 1, we can see that the largest relaxations are for the metal atoms in the surface layer, 
-5.7% and -3.4% for the Sr-terminated and Ti-terminated surfaces respectively. The outward 
relaxation of the second-layer Sr atom on the Ti-terminated surface is also noteworthy. The 
surface-layer oxygen atoms show almost no displacement on the Sr-terminated surface. 

Previous surface structural refinements have been carried out by Bickel et al. [21] and by Hikita, 
Hanada and Kudo [20] using (LEED) and (RHEED), respectively. In Table 2, we compare our 
structural parameters with the two experimental ones. Both experimental groups assumed that 
the oxygen and metal atoms remain coplanar in the second and third layers in order to simplify 
the refinement procedure; since the scattering strength of O is much smaller than that of Sr 
and Ti [21], our theoretical comparison is made to the position of the subsurface metal layers. 
Thus, we define A<i 12 as the change (relative to bulk) of the first interlayer spacing, as measured 
from the surface to the subsurface metal z coordinate, and similarly for Ad 23 . The quantity 
s measures the outward displacement of the surface-layer oxygens relative to the surface-layer 
metal atoms. 

It is evident that the agreement between the theory and the experimental refinements is not 
very good. This is not surprising in view of the fact that the experimental structures are in 
poor agreement with each other [30]. Our values for the rumpling s of the surface layer are in 
rough qualitative agreement with the experimental ones, though somewhat larger in magnitude. 
However, we predict a substantial contraction of the interlayer spacing cq 2 , while in most cases 
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the experimental refinement indicates an expansion instead. The one exception is for the Sr- 
terminated surface, where our theory is in reasonable agreement with Ref. [21], although not 
with Ref. [20]. The agreement for the second interlayer spacing can be seen to be mixed. 

The disagreement between the two sets of experimental numbers [30] suggests that the experi- 
mental refinements should perhaps not be taken too seriously. The expected quality of a LEED 
or RHEED refinement is not well established for a complicated metal oxide surface such as this 
one. In the work of Bickel et al., [21] the authors were not able to determine independently the 
proportions of the surface exhibiting the SrO and Ti0 2 terminations; they assumed that the 
two appear in equal proportions. Refining the structural parameters for both simultaneously, 
they then obtained an /2-factor of 0.53. While this was argued to be "acceptable in view of the 
complexity of the structure," it nevertheless seems uncomfortably large. In the work of Hikita 
et al., [20] the surface was prepared in different conditions in order to obtain SrO and TiC>2 
terminations separately; for these, i?-factors of 0.28 and 0.26 were obtained, respectively. While 
this would thus appear to be the more reliable experiment, unfortunately it is in no better 
agreement with the theory than the refinement of Ref. [21]. One possible problem with both 
experimental refinements could be the arbitrary assumption that there is no buckling in the 
second metal-oxygen plane. Our results indicate that there is a substantial buckling in that 
layer, especially in the case of the subsurface SrO layer on the Ti02-terminated surface. 

Also included in Table 2 are theoretical estimates of Mackrodt [29] using an interatomic po- 
tential based on Kim-Gordon pair potentials. This theory gives results that are qualitatively 
similar to ours, although the pair-potential model appears to overestimate the size of the inter- 
layer relaxations and underestimate the degree of surface rumpling. The pair-potential model 
also predicts a significant rumpling in the subsurface SrO layer of the Ti02-terminated surface, 
similar to what we reported above in Table 1. 



4 Testing for a ferroelectric monolayer 

Bulk SrTiOa is an incipient ferroelectric: it nearly becomes ferroelectric at very low temperature, 
and is apparently prevented from doing so only by quantum zero-point fluctuations [31]. In our 
work on surfaces of FE BaTiOs, we found some tendency for an enhancement of the FE order in 
the first few surface layers, especially for the case of the Ti02 termination. Thus, it is intriguing 
to speculate that if a similar tendency exists in SrTi0 3 , it might lead to the formation of a FE 
surface layer on top of a paraelectric bulk material, at least at low temperatures. In fact, the 
presence of just this kind of FE monolayer at the surface has been predicted by Ravikumar, 
Wolf and Dravid [18] for the SrO termination of SrTi0 3 . However, this prediction is based on 
an empirical interatomic potential developed for ionic systems, and it is not clear how far such 
a pair-potential approach can be trusted for deciding such a delicate question as the appearance 
of surface FE order. Thus, we have undertaken to check whether such a FE surface phase might 
occur in the context of our ab-initio calculations. 

Note that while Bickel et al. have discussed the rumpling of the surface layer in terms of 
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a "ferroelectric relaxation" normal to the surface [21], such a rumpling does not break any 
symmetry of the ideal surface, and is not qualitatively different from the rumpling observed on 
other oxide surfaces. Thus, when we speak of a FE surface layer, we shall restrict ourselves here 
to symmetry-breaking distortions, i.e., frozen-in displacements parallel to the surface. Indeed, 
Ravikumar et al. predict enormous displacements of this kind for the SrO-terminated surface: 
Sr and O shift by 0.43A and -0.33A, respectively, along a (100) direction, with a resulting 
reduction in surface energy of 0.11 eV/cell. When the primitive lxl periodicity of the surface 
cell is allowed to be broken, they find an even lower-energy c2x2 structure. 

With this motivation, we have investigated carefully the stability of our lxl surface with respect 
to displacements of the surface atoms along the x axis, using the reduced set of symmetries 
as discussed in Section 2. Starting from the structure of Table 1, an additional FE distortion 
along x is imposed, the value of this distortion being taken the same as for bulk BaTiC>3 as 
computed in Ref. [24]. Then we allow this structure to relax fully. We did these simulations for 
both kinds of termination, and for both theoretical and experimental lattice constants. 

For the case of experimental lattice constant, the resulting relaxed displacements along the 
x direction, in units of lattice constant, are as follows. For the SrO-terminated surface, the 
first- and third-layer Sr displacements are 0.005 and 0.006; the second- and fourth-layer Ti 
displacements are 0.005 and 0.008; and the oxygen displacements are -0.007, -0.002, -0.006, and 
-0.004 in layers 1 through 4 respectively. For the TiCVterminated surface, the first- and third- 
layer Ti displacements are 0.017 and 0.016; the second- and fourth-layer Sr displacements are 
0.010 and 0.011; and the oxygen displacements are -0.011, -0.008, -0.007, and -0.009 in layers 1 
through 4 respectively. The total energies lie below those of the high-symmetry structure (i.e., 
no x distortion) by 0.001 and 0.026 eV for the SrO- and TiCVterminated surfaces, respectively. 
We see that the final FE distortions left after the relaxations are indeed very small, and are not 
much greater in the surface layer than in the deeper layers. When using the theoretical lattice 
constant instead of the experimental one, we found only a very weak distortion of energy 0.001 
eV for the Ti0 2 -terminated case, and no observable distortion for the SrO-terminated case. 

These results indicate that there does exist the possibility of a small FE surface distortion 
for both surface terminations of SrTiOs. However, the double- well depth that we found for 
the SrO surface is very much smaller than the one predicted by Ravikumar et al. [18]. We 
actually find a somewhat larger tendency for FE distortion for the Ti0 2 surface, but still the 
distortion amplitudes (maximum 0.07A) and energy (0.013 eV/surface) are nearly an order of 
magnitude smaller than the SrO-surface prediction of Ref. [18]. Our predicted FE distortion is 
small enough that it may easily be destroyed by thermal (or even quantum) fluctuations, and 
it seems unlikely that it would appear at room temperature. However, it seems possible that 
the FE distortion might survive and be observable at very low temperatures. 
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Fig. 1. Grand thermodynamic potential F as a function of the chemical potential n<rio 2 , for the 
two types of surfaces. Dashed and solid lines correspond to type-I (SrO-terminated) and type-II 
(Ti02-terminated) surfaces respectively. 

5 Surface energies 

We discuss next the surface energetics. To compare energies of surfaces having different stoi- 
chiometry, we follow the approach described in Ref. [13]. That is, we consider each surface slab 
to be built from TiC>2 and SrO units, and compute the grand thermodynamic potential F for 
each type of surface as a function of the Ti0 2 chemical potential ^Tio 2 - The energies of the bulk 
crystals of SrO and Ti0 2 were calculated with the LDA using the same pseudopotentials, and 
the same pane- wave cut-off. The final results are shown in Fig. 1. It can be seen that both sur- 
faces have a comparable range of thermodynamic stability, indicating that either termination 
could be formed depending on whether growth occurs in Sr-rich or Ti-rich conditions, with a 
very small preference for the TiC>2 termination. 

The average of F for the two types of surface, which we shall denote as i? sur f, is independent of 
//ticv Thus, this quantity is convenient for comparisons. The value we found for E SUT { for the 
SrTiOa (001) surfaces is 1.26 eV per surface unit cell (1358 erg/cm 2 ). This is very similar to 
the corresponding value 1.24 eV that was obtained for the BaTi0 3 (001) surfaces (in the cubic 
phase). 

To compute the surface relaxation energy -E re iax, we computed the average surface energy -E^m-ei 
for the unrelaxed slabs (i.e., atoms in the ideal cubic configuration), using the same k-point 
sampling as for the relaxed systems. We obtained -Eunrei = 1-44 eV. Thus, the relaxations 
account for 0.18 eV of the surface energy per surface unit cell, accounting for around 15% of 
the total surface energy. 
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r x M r 

Fig. 2. Calculated electronic band structures for SrTiC-3. (a) Surface-projected bulk band structure, (b) 
SrO-terminated surface (slab I), (c) TiC>2-terminated surface (slab II). The zero of energy corresponds 
to the bulk valence band maximum. Only the lowest few conduction bands are shown. 

6 Surface electronic structure 

We focus now on the LDA calculated electronic structure for the surfaces slabs. Although the 
LDA is well known to underestimate the band gaps, we can be assured that the results given 
here are at least qualitatively correct. The band gap we obtained for bulk is SrTi0 3 1.85 eV, 
to be compared with the experimental value of 3.30 eV; this level of disagreement is typical for 
the LDA. 

Not surprisingly, our results for the band structures of the surfaces of SrTiOs are very similar 
to the case of BaTi0 3 surfaces [13]. The computed LDA electronic energy band structure is 
given in Fig. 2 and the direct energy gaps are given in Table 3. We see that the band gap for 
the SrO surface almost does not change with respect to the bulk value, and no in-gap state 
occurs. For the Ti02 surface, there is a substantial reduction of the band gap. However, there 
are no deep-gap surface states, in accord with experimental reports [6]. 

From Fig. 2(c), we see (as is also true for BaTiC^) [13] that there is a tendency for valence-band 
states to intrude upwards into the lower part of the band gap for this surface, especially near the 
M point of the surface Brillouin Zone, that is the top of the valence band. In order to analyze 
the character of the valence-band state at the M point of the Ti0 2 surface, the charge density 
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Table 3 

Calculated electronic energy band gaps for the relaxed surface slabs (eV). 



Slab 


Band gap 


Slab I 


1.86 


Slab II 


1.13 


Bulk 


1.85 



at this state was obtained. The results are very similar to those shown in Fig. 4 of Ref. [13] for 
the case of BaTiC>3. These states correspond to 2p orbitals lying in the surface plane and 
having very little hybridization with Ti 3d levels because of the existence of four nodal planes 
[(100), (110), (010), and (110)] intersecting at the Ti sites. 



7 Summary 

In summary, we have carried out LDA density-functional calculations of SrO- and TiCv 
terminated (001) surfaces of SrTiOs. By minimizing the Hellmann-Feynman forces on the atoms, 
we obtained the optimal ionic positions. Previous experimental LEED [21] and RHEED [20] 
surface structure determinations are found to be significantly at variance with our predictions, 
as well as with each other. 

Our calculations do not support the existence of a large FE distortion in the surface layer of the 
Sr-terminated surface, as had been previously proposed theoretically [18]. While we cannot rule 
out a weak FE surface instability, especially for the Ti-terminated surface, our results indicate 
that it would be sufficiently weak so as to be easily destroyed by thermal fluctuations except 
at quite low temperatures. The equilibrium state of the surface was obtained as a function of 
the relative Sr and Ti stoichiometry. The average energy of the Sr- and Ti-terminated surfaces 
was found to be ~ 1360 erg/cm 2 , and the surface relaxation energy was calculated to be ~ 190 
erg/cm 2 . In agreement with experiments, no deep-gap surface level is found. However, as found 
for BaTiC>3, there is a significant reduction of the electronic band gap of the Ti02 surface. 
This reduction arises from surface states in the lower part of the bulk band gap which can be 
regarded as resulting from the intrusion of valence-band states of O 2p character, particularly 
in the vicinity of the M point of the surface Brillouin zone. 
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